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Abstract 

We present explicit analytical formulae for two-body matrix elements between Pauli-projected 
single-particle orbits generated by gaussians shifted from the centre of a nuclear core, which is 
populated by nucleons occupying harmonic-oscillator orbits. Such shifted gaussian orbits appear in 
the generator-coordinate model of two-cluster systems and, as ingredients of generating functions, 
in the resonating group model as well. A Pauli projection is to be understood as a projection 
onto the subspace not occupied by the core nucleons. This is an exact procedure to take into 
account the Pauli principle in the limit in which the core is infinitely heavy, thus the procedure to 
be presented is very useful in the microscopic description of large-core plus small-cluster systems, 
such as involved in a-decay problems. 



31. Introduction 



The more complex the nuclear system, the more complicated it is to allow for the Pauli principle, 
single-configuration models, like the extreme shell model, the Pauli principle is exactly allowed for 
■py choosing the single-particle (s.p.) orbits orthogonal to each other and prescribing single population; 
nence the interpretation of the Pauli principle as an exclusion principle. For multiconfiguration one- 
centre models with a passive and massive core, like the multiconfiguration shell model, the Pauli 
principle between the core and valence particles is exactly observed by choosing the valence orbits 
orthogonal to those occupied by the core particles. However, for a multicentre system, such as a 
two-cluster system, the antisymmetrization between the clusters has to be carried through explicitly, 
and that implies tremendous technical complications [||. 

In this paper we consider the case when one of the clusters, the "core", is substantially heavier 
than the other, "the cluster". In the limit, in which the core is completely passive and infinitely heavy, 
the Pauli principle can again be taken into account by an orthogonality requirement, but it is not so 
easy to fulfil this requirement. The subject of this paper is just to show a technique for the treatment 
of the most critical aspect of such a model, that of the calculation of the two-body potential matrix 
elements. 

The wave function of the core A is taken to be an harmonic-oscillator (h.o.) shell model Slater 
determinant, and that of the cluster a, which is assumed to contain at most four nucleons, is restricted 
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to be Os h.o. orbits around a centre displaced from the centre of the core. The two-centre shell model 
underlying this picture forms a basis of the generator-coordinate model (GCM) of the two-cluster 
system. The orthogonality condition is imposed on the s.p. functions of cluster a by an operator that 
projects out the states filled in the core. 

Though this model may look rather restrictive, it can serve as a starting point for more sophisti- 
cated models. In particular, the condition of the infinite mass of the core can be lifted by applying a 
transformation that shifts the centre of mass (cm.) of the system from the centre of the core to the 
joint cm. By further transformations the GCM can be replaced by a resonating-group model 
or by one in which the relative motion is represented by an h.o. basis |^ or by a single-centre basis 
spanned by gaussians of different widths [Q] . Any of these transformations only involves the generator 
coordinates (and not the physical ones), so they do not interfere with antisymmetrization and can be 
performed on the matrix elements. All of these transformations can be performed analytically, but 
they should be performed before the angular-momentum projection. 

The model with Pauli-projected s.p. orbits is useful especially for the description of core-|-a 
systems. Since the relation of the Pauli-projected orbits to the core orbits is the same as that of the 
valence shell-model orbits, the Pauli-projected orbits can be treated, to some extent, as shell-model 
orbits. In particular, the configurations they can form are just like ordinary shell-model configurations. 
The only difference is that they are not orthogonal to the ordinary valence orbits, thus a diagonalization 
over a non-orthogonal basis is required. Such cluster configurations have been used in the description 
of Q decay |, |, §. 

To calculate the overlaps and s.p. matrix elements between such Pauli-projected s.p. orbits is 
straightforward with the well-known techniques Q. The calculation of the two-body matrix elements 
is, however, not so simple. The two-body matrix elements reduce to integrals of products of h.o. 
and shifted gaussian functions. It seems to be natural to expand the gaussians in terms of h.o. 
functions Q and then use the conventional technique, but this recipe leads to infinite summations, 
and the numerical expenditure is enormous. The procedure we show avoids this, by evaluating these 
integrals directly. 

The plan of the paper is as follows. In §2 we briefly summarize the essentials of the method of 
the Pauli-projected s.p. orbits. In §3 we express the matrix element of an interaction term as a sum 
of sixteen terms and derive a general formula for these terms. In §4 we specialize this formula for 
some terms, and in §5 we give a brief summary. The paper is concluded with an Appendix containing 
a glossary of the formulae we used in the derivation. 



2. Pauli-projected single-particle orbits 



A general state of the a + A two-cluster system can be written in the form ^ = A{(f>LM{T^) 
^ai(,a)^A{£.A)}, where ^ are antisymmetrized cluster internal states, the function describes the 
relative motion of angular momentum LM and A is an intercluster antisymmetrizer normalized in the 
pattern of ^ = A^-^s J2iLi{-'i-)^'Pi- A GCM representation can be obtained by expanding (Plm{t^) 
in terms of angular-momentum projected gaussians displaced with respect to the centre: 

4>LM{r) =J2fk(iSkYLM{h)Vs„{T^), (1) 
k •' 
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where the gaussians are Os h.o. functions of parameter b, 



V^s(r)=fzV^'exp(-^(r-s)n. (2) 



Accordingly, ^ will take the form ^ = J2k fk^k, with to be discussed below. 

Let us assume that the two intrinsic states are single configurations of h.o. states ^qqq and 
'PniUmi^ (ni, /i, mi, n^, "t-a), of parameters a and /5, respectively. If, in addition, the core is 
taken infinitely heavy, <I>a reduces to a Slater determinant of Wnikmi} choosing b = 4a, so 

does <fg^{r)^a{S,a)- The corresponding a + A-particle state thus reads 

= ^{y^'"(r-Sfc)$„(e„)$e(ec)} (3) 
= [(a + A)!]-i/2det{<(ri)xi(l)...<(r„)xa(a) 

X^lhmi (ra+l)Xa+l (a + {ra+A)Xa+A{a + A)}. (4) 

The coordinates rj are those with respect to the total cm., but, the mass being infinity, the cm. can 
be fixed at the origin, so may be interpreted as the usual s.p. coordinates. The symbol Xi(i) stands 
for a spin-isospin state. The function may be cast into an even more useful form, 

= [(a + ^)!]-V2det{V;s,(ri)xi(l)...V's,(ra)Xa(a) 

hmi i^a+l 



where 



Mr,) = (1 - VjMir,) = ^^(r,) - ^ ^f,,^,(r,)(<,^J^^), (6) 



i=l 

with Vj = J2f ^niiimiirj){^niiimi\- '^^^^ ^-P- statcs V's(r) Pauli-projectcd s.p. orbits. The 
function ^'^ of eq. @ is equal to that of eq. @ because (||) is obtained from (Q) by adding linear 
combinations of columns a+l,...,a + ^to columns 1, a. 

The advantage of using (|5|) rather than (^) is that the s.p. states involved in eq. (H) are all 
orthogonal to each other. The matrix elements between two states of the type of eq. (|5|) can be 
expressed in terms of the s.p. states more easily, which helps to separate the problem of the valence 
nucleons from that of the core [^. In particular, the two-body matrix elements of Slater determinants 
built up of orthogonal s.p. states can be expressed simply by the two-body matrix elements of the 
s.p. states. 

S3. Matrix elements 



We shall consider the two-body matrix elements of the form factor, e 2''(''2 '"i) ^ of a gaussian 
two-body potential between the Pauli-projected s.p. functions V's(r): 

(V's(ri)V's(r2) |e-^'^('-^-^i)'| Vs'(ri)V's'(r2)) 
= ((1 - Pi) V9^(ri) (1 - V2) ^^ir2) le-H"-^-^^)'! (1 - V,) ^-,(ri) (1 " ^2) ^^'(rs))- (7) 

It should be noted that the s.p. states in the bra and ket have no definite angular momenta, but a 
state of good angular momentum can be easily constructed from any of V's by taking /dsy;m(s)V's- 
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Since, however, the projection of the total angular momentum in eq. (||) also involves the s.p. states 
other than those in the two-particle matrix element (|^, it is more natural to include this projection 
at a later stage, which is out of the scope of this paper. Moreover, as is obvious from eq. (|I|), the 
projection of the full orbital angular momentum only requires a single operation by JdsYcMi^)- The 
only thing we shall do in this paper to facilitate angular-momentum projecting is that we shall express 
the matrix element in terms of s, s, s' and s'. 

By letting 1 — Vi act term by term, this matrix element is written as a sum of 16 terms. Each of 
these terms contains two-body matrix elements of different combinations of h.o. and gaussian wave 
functions. To write the two-body matrix elements of the orbits ip'^ and V^p^ as special cases of a 
general formula, we introduce the notation 

M{qiq2\qW2) = (^'Vs (rO^^IVs (rs) je-^'^^^^-r^)' | <,(ri)Pf ^^(rs)), (8) 

which runs over all terms when each quotient Qi runs over and 1. With these, (|) can be expressed 
as 

{i;s{ri)Mr2) |e"5''(r2-ri)^| V;s'(ri)V's'(r2)) = M(00|00) - M(10|00) - M(01|00) 
-M(00|10) - M(00|01) + M(11|00) + M(00|11) + M(10|10) + M(01|01) + M(10|01) 
-FM(01|10) - M(01|ll) - M(10|ll) - M(ll|01) - M(ll|10) + M(ll|ll). (9) 

It is obvious that any two-body matrix elements that involve shell-model orbits , ^ instead of some 
of the Pauli-projected excentric states ips are built up from the same building blocks, but are simpler, 
and can be derived analogously. 

The calculations are simplified by the following symmetry relations: 

M(10|ll) = M(01|ll), M(ll|10) = M(ll|01) = M(01|ll)', (10) 

M(00|10) = M(00|01), M(10|00) = M(01|00) = M(00|01)', (11) 
M(10|10) = M(01|01)', M(11|00) = M(00|11)', M(10|01) = M(01|10)', (12) 
where M{qiq2\q[q2y differs from M(gig2|(?i(?2)' ™ that the parameters s and s' are exchanged. 
It is convenient to factorize the h.o. function as 

With the introduction of Ani{s) through 

{^tWllm) = Anl{s)YU^), (14) 

the matrix elements of eq. ^ take the form 

M{qiq2\q'A) = (sWg^ {sWg[ {s')M,^, {s) 

X ^"I'l i^)^hnii (8)^2/2 {s)Yl^m2 (s')>l;m', {^'T K'^V., {s')Yv^m'^ (s')* 

{n'A'.m'.} t 
'-II l' q. 

I 

xjdr^j rfr2<,^„^(r0*^^,,,„,(r2)*e-^i:li(-^(r'-)^+'3»^^) 

,,-^.(r.-r.)^^J (r^)^^ (r^)^-|EL.("^(r.-s')^+/^^.n, (15) 
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with the notations = (1 — qi)a, Pi = QiP and 



a 



X 3/4 



1-9 



(16) 



As to the summation hmits, if = 1 then {nj/jmj} runs over the filled orbits, but, if qi = 0, then 
there is just one term, {rii = 0, /j = 0, rrii = 0}. For Ani, an explicit formula is given in eq. (53) of the 
Appendix. 



It is useful to couple the angular momenta of the h.o. wave functions of the same arguments and 
make use of the formula [^] 

^ ^ n 

(the coefficient T^[i-^^n2i2 ^® given in eq. ( |5^ ) of the Appendix), to obtain 

//3\3/2 

M{qiq2\qW2) = [ -) ^,.{s)M,,{s)M,,^{s')Mq>^{s') ^ A„,j,(s)A„,i,(s)^„,;, (s')^„,;, (s') 



{n'i'} , 



AA'LAf AfAf' 



with 



jNN' 



(18) 



dri / dr2 



LM 



X e~5Eli("»(r»--)'+ft'-?)e-5'?(r2-ri)%-fELi(":.(r.-«')V/3H). (19) 

The summations over the angular momenta are restricted by the well-known rules, i.e. \li — I'll < X < 
h + I'l, and M = —L, L. The coefficient T^[i^n2i2 non-zero if a triangle inequality (see (|55|) in the 
Appendix) of the quantum numbers is satisfied. In eq. (18) the spherical harmonics can be recoupled 
so that those belonging to the same argument can be combined: 



y;,(s)y,,(s')] k(s)i^^(s') 



LM 



h I'l A 



I I' L 



LAI 



(20) 



where J = \/2J + 1 and 



(/1O/2OI/O). 



(21) 



In §4 we shall see that at least one of the angular momenta in the 9j-symbol is zero in all practical 
cases, therefore ( pO[ ) reduces to a simpler form. 



The integral IIIx'')lm determined using the the generating function of ^.^ P , 



(22) 



With this, (f^i^ can be expressed as 



3/4 



Cnl 



^2n+l 



IT J (2n + 0!dp2'^+' 

Cnl = (-1)' 



n 



dpYim{p)g{p) 

p=0 

!(2n + 2/ + l)!!y/2 



47r2"+' 



(23) 
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Thus the integral I^x\')lm determined from the matrix element of the generator function 

G(p,p') = j dv^j dr2g{jy)e''2T.lA-ri-^--?+Prrl)^-\^^^^^^ (24) 



as 



tNN' _ ( Cn\CN'\ 



(AA')LM y^j (2N + \)\{2N' + \')\ 

' ' dpj dp' [Yx^{p)Yy^,{p')]^j^,jG{p,p': 



Qp2N+\Qpl2N'+\> 



p=0 
p'=0 



(25) 



With the abbreviations 



d = iaia2-^^ Ai = J, M = ^, P = 2V^p, P' = 2^p' , 

Bi = - Ai, B2 = j0-A2, B = 2, V = 2Aiaia[ + 2A2a2a2 + B{aia'2 + a'ia2), 
h{ai + Q2) — AiOi — ^2«2 ~ Baia2, u' = \{oi'i + a'2) — Aia'-^ — A20i2 ~ Ba'ia2, 



u 



(26) 



t = {2Aiai + Ba2), r = {2Aia[ + Ba^), t' = {2A2a2 + Bai), r' = {2A2a2 + Ba[), 
Q = ts + rs', Q' = t's + r's', 

the double integral G{p, p') is expressible as 

/a 2 \ "^/^ 

G{p, p') = expl-us^ - u's'^ + vss' - BiP^ - ^2^'^ + ^PP' + PQ + P'Q'}. (27) 



The explicit forms of these coefficients in the exponent for different values of qi and q'^ are collected 
in table 1. To be able to carry out the operations prescribed by (p5|), we expand the exponentials 
containing P and P' into the power series 



G(p,p') = exp{-ns2-nV' + W} V [ ^} , , „2.i+^i+y 2.2+^2+^ 

V ^ / .i.2m,..K^i!^2!w!/^2!K! 

A1A2/ 

X ni™i(p)*nirn,(Q)n2m2(p')*n2rn2(Q')^Zrn(p)i1™(p')*, (28) 

mim2m 

where we used the expression 

^ ^ A.-TTT] f 

(ab)" = a-b^Y^Cni E YUmmihT, C^i = _ ^),,(^ + / + i)n ' (^9) 

In this formula J2i{n) stands for a summation over I = 0,2, ...,n — 2, n if n is even, and I = 1,3, ...,n — 2, n 
if n is odd. 

It is straightforward to recouple the spherical harmonics in (p^): 

E Yx,m, {vYYx.m, {Q)Yx,m, {vYYx.m, (Q')^m (P)>;m (p')* 

m\m2m 

= E {-l)'y^W{\^Ll^;\2\)SliSli [Yx{p)Yx>miM [n,(Q)n,(Q')]^^, • (30) 

AA'LAf 
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After substituting (30) into ( [2^ ) and inserting ( |2^ ) in ( [25| ) the derivations and the integrations in (|2 
can be carried out, with the result 



2Ar+A27V'+A' 



Ml=0 /i2=0 Ai(/ii) A2(/i2) 



(31) 



where 



{NXN'X')L 



2N+X+2N'+\' 



Cat A CAT' A' C^i Ai (^^12 A2 



X E / , , , ijr-B2^-B"E(-l)'^(^i-^^^^^2A)C,,g,\,5,\„ 



and the summation indices should fulfil the constraints 

2i^i+ Hi + K = 2N + A, 2i/2 + Ai2 + K = 27V' + A'. 



(32) 



(33) 



The vectors Q and Q' are linear combinations of s and s'. To facilitate the angular-momentum 
projection, it is useful to expand the Q- and Q'-dependent factor of eq. (|3l| ) in terms of s and s': 



Q^'Q'^' [n,(Q)>A2(Q') 



LM 



fJ-l M2 



E E t^V^'i^'^it"^V^2-'=2gfci+fc2g/m+^^2-fci-fc2 ^ F^'J'^^i^^f^ [Yx'^{s)Yx'^{s')\ , (34) 



fci=0fc2=0 



A'^Aj 



where the coefficients F'x^x^x^x^j^ are given in eq. ( |57| ) of the Appendix. One can verify this expansion 
immediately by using the well-known expression Q 



(a + bTYU^ + b) = E E E ^AfA2/ [^i (^)n2 (b)' 

fe=0 Ai(A:) A2(n-fc) 



Im 



(35) 



where H'x^x2i given in eq. ( p6| ) of the Appendix. 

With ([3T|) and (p4|) substituted into eq. ([T8|), our final expression for M {qiq2\q'iq'2) becomes 



M{qiq2\qiq2) = [^^j Nq,{s)Nq,{s)^^q[{s)^^^'^{s) exp{-us2 - u's" + W} 

/ /i A 

X E ^n,Zi(s)^n2;2(^)^n;Zi(s')^n^«^(s')EEE(-l)''^'^"^"'^^'| ^2 A' \ Sli,S{, 



XX' W L 



I I' L 



I l' q. 



2N+X2N'+X' 

ErpNX rpN'X' MNXN'X')L 



Ml M2 



Mi=0 M2=0 Ai(/xi) A2(m2) 



^ ^ ^ ^ki^tj,i-ki^lk2^ffi2-k2 gki+k2 g'fJ.i+fJ.2~-ki-k2 
fci A:2 



E ^a7a'2 w E(-i)^'^(^^'^'i^2; ^^O^^A^^^;^ E ^^^^^^^ 

A' A' L' M' 



(36) 
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§4. Special cases 



In this section we substitute the values of qi and q'^ into eq. ( |3q ) and determine the concrete 
forms of M {qiq2\q[q2) • 

(1) M(ll|ll). This case reduces to combinations of two-body matrix elements of h.o. functions. 
These matrix elements are extensively used in shell- model calculations 0|. From table 1 we see 
that t = t' = r = r' = 0, i.e. Q = 0, Q' = , n = n' = f = 0, and then only the terms with 



^J'i = ki = \ = 0, {i = 1,2) survive in equation (36). Moreover, as 



nOOOO 



^'\i\2\[X'2L = '^AiO'^AaO'^A'iO'^A^O'^LO (37) 

the the Racah, the 9j-symbol and the S^, coefficients in equation (^) are reduced as: 

W{im;OL') = ^ L , (38) 

h 1'2 A' =6xx'6u ^ ' W{hl[hl'2-Al) (39) 

I I' 



and 



By using these results, eq. ( [3q ) reads as 

M(ll|ll) = (n ^n,h(^)^n./.(5)A„;«;(s')A4z^(s')E^(-l)V^(^i^i^2Z^;A0 



^1 ^ .^I'lV-/ n^i^x- , iL^i^y- , ^ -y , ■ 

{n'j'} 

m NN' 



In this case the coefficient Gq^q^ ^^'^ can be expressed in a simpler form. From eq. (|3^ ) it follows 
that K = 2N + A — 2z^i and z^2 = N' — N + vi, and then by substituting Bi, B2 and B (see table 1) 

into eq. (|3^), reducing the Racah and the 5;^ coefficients, Gq^q^ ^'"^ is given by 

^ (4.)^ {^f c„.c„.. (2^,^^^^^"'^^ (-1)— A (^)""'" 

2N+\-ui 

^ ui\{N' -N + vi)\{N - vi)\{2N + 2A - 2ui + 1)!! ' ^^^^ 



^1 



The summation over vi can be carried out Q and our final expression is 

Gqooo = j CNXCN'X J (-1) 2 (47r) A 

(2Af + 2A^' + 2A + 1)!! 



Af!Af'!(2iV + 2A + iy.\{2N' + 2A + 1)!! ' 



(42) 
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(2) M(01|ll). Now we do not have so simple expressions as in the previous case. The first simphfication 
is that r = r' = 0, i.e. Q = is, Q' = t's and u' = v = (see table 1). The summation over ki{i = 1, 2) 
is then restricted to ki = = 1,2). The coefficient FxiX2X'^x'^L easily be evaluated: 

KtATA^L = ^^KL^X'^A. ■ (43) 

From this it follows that A2 = 0, thus the Racah coefficient and S^^, in equation (36) can be written 
as: 

Will'LO; LL') = = -^6vl' (44) 

An other simplification to be used is that when q\ = 0, then ni = li = 0. From this it follows 



that the coefficient Tj^, „, is 



NX 

'I'l 

NX 



T00n[l[ = ^n[N\x (45) 

and the 9j-symbol, similarly as in the first case, can be expressed by a Racah coefficient. 
Let us put all these together to find 

M(01|ll)= f^Y\^p{-us'} ^ A^^i^{s)A^,i,y)A^,i,{s') J2 X'L{-iy'W{L%LX';l[h) 

{n'!'} 

m' N' MiAt2AiA2 

(3) M(01|10). As we can see in Table 1, all coefficients differ from zero. There is one gaussian in 
the bra and one in the ket, that is qi = q'2 = and ni = li = n2 = I2 = 0. It follows that 



'^00n[l[ - ^n'^N^l[Xi ^^2^200 " ^n2A"^«2A', (46) 

and the 9j-symbol, as it contains two zero elements, can be substituted by 

^ A \ (_iyi'+i2-L 

h A' = ^-^2 ^V^x^ih^i'V^h'h (47) 

\l V L ) Ijl'i 



That is all one can do in this case and thus equation (^6D becomes 



/ x3/2 2n[+l[2n2+l2 

M(01|10)= (-] exp{-ns2-nWW}E^n2Z2(s)^n'/,(5')E E E 

"2*2 L Ml=0 M2=0 



•I'l 



E\ ' /-("itin2(2 \ " \ " ^fei^^ti-fci ./fc2^'M2-fc2gfcl+fc2g'W+A'2-fcl-fc2 \ " p^llkl^l2k2 
/ ^ M1A1/12A2 / ^ / ^ / ^ A1A2AJA2-L 

Ai(/xi) A2(m2) A:i=0fc2=0 X{X'2 

L2E(-l)'^^''^^'w^(^2/iA;A'2;LL05,^;,5,f,,En'AfKs)n'M'(sO*. 

L' M' 

(4) M(01|01). Now qi = q[ = {) and m = h = n[ = l[ = 0, consequently T^^q = 5nx and 
therefore = ki = 0. Then the coefficients will reduce to 

^00^42^2 _ rrM2fc2 //^o^ 

^OA^AIA^L - 7|^'^^A2^A;Ai^A2' (48) 
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and after repeating the steps used in the previous cases, i.e. reducing the 9j-, Racah and S-symbols, 
the expression (|36|) yields 



/ 2/o\3/4 2N'+L 

M(01|01)= (^) eM-ns''-u's'^ + vss'}Y.An,iMMi'S^')Y.T.<Ui'. E 

L N' fj.2=0 



n' l' 



xSt^,,Sl/,,Y.^L'M'{s)YL'M'{sr, 



M' 

where the coefficient G can be written as 

^(OO^'A')A' _ /4/37r\3/2 , .2^'+A' (-1)-^ 

G^oom2A2 -l^j CN'x'[2VP) 2. ——B^C^.y, (49) 

(5) M(00|11). In this case both gaussians are in the bra, and u' = v = 0, r = r' = 0, therefore 
Q = ts, Q = t's just as in case (2). In addition, now ni = li = n2 = h = 0, so the result of case (2) 
can be simplified further: 



/ .3/2 27V+A 27V'+A' '/.I /M2 ,,,,,,,,, 

^(00111)= (^) exp{-ns^}EA„;.i(.)A„,,,(.')E E E E E ^l^felf 



n'^l[ L /il=0 /i2=0 Al A2 

"2*2 



>,iMii'M2^m+/.2 5L ^ YLM{s)YLM{s'r 

M 

(6) M(00|01). Using the results of case (4), one finds 

^3 \ 3/4 ^ 2^+(2 

^2 ^ 1 nA' c^'2 



00^2 «2 

rijij ^^2=0 k2 



M(00|01)= (^) exp{-ns2 - ^xV' + W} ^ £ ^(00njp'^^i/fc2/M2-/c2^/c2^/M2-fc2 

i, 

A', Ao 1 '"2 



E ^rfA\\>T(-l)'^4\; E>^A>.(s)rA>.(s')*. 



'1''2 



3/2 

M(00|00) = ( ] exp{-us2 - u's'"^ + vss'} (50) 



(7) M(00|00) This is the simplest case, overlap of four gaussians is just 

a 

§5. Conclusions 

All summation in formula (p^) are finite, so after angular momentum projection it is straightfor- 



ward to apply it, or its special cases given in §4., in numerical calculations |1C]. If, in any of the states 
involved there are other nucleons in Pauli-projected orbits, their overlaps may bring in factors of the 
form of exp{u;ss'}. Grouping these factors together with exp{t;ss'} of eq. the product has to be 
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expanded into its multipoles. The functions Y\m{s)Y/<^m{s')* appearing in this way should be coupled 
with yL'M'(^)^'M'('S')*j respectively. The infinite summation involved in the multipole expansion 
is then eliminated by the angular momentum projection, which can be achieved by the operations 
JdsYljy^{s) and J ds'Yc'M' is')- 

I am indebted to Prof. R.G. Lovas for carefull reading of the manuscript. 

This work was supported by the OTKA grant (No. 3010). 



§6. Appendix 

In this paper we used the standard definition of the harmonic oscillator function (see e.g. Q) 

with L^^^^^(x) being associated the Laguerre polynomial 

_ (2n + 2/ + l)!! " {-ir2^ 
" 2- (2s + 2/ + l)!!s!(n-s)!^ ■ ^ 

The radial part (see eq. (|l^) of the overlap of a harmonic oscillator function of width parameter 
/? and a gaussian of width parameter a and of displacement vector s is 

/n!(2n + 2/ + l)!!\^/V 4a/3 X^/^ 

fr'Q ^ k\{n-ky.{2k + 2l + l)\\\a + f3j [a + (3 J ^ ^ 

The combination coefficients of product of two h. o. functions of the same argument ae given by 

i 

tNL _ ( n\n'l2^^^{2n + 21 + l)!!(2n- + 21' + 1)!! ^/^ W , 

^nin'v - y (2iV + 2L + l)!!2"+'2«'+''iV! ) "l^")^ 

^ {2s + 2s' + 1 + 1' + L + l)\\{s + s' + \{l + l' - L))\ 

^ (2s + 21 + l)!!(2s' + 21' + l)!!(n - s)\s\{n' - s')\s'\{s + s' - N + \{l + I' - L))\ ' ^ ' 

where the quantum numbers must satify the triangular condition 

\2n + l-2n' -1\<2N + L<2n + l + 2n' + 1'. (55) 



The expansion coefficients in equation ( |35D are 



Tjnk _ I Ck\^Cn~kX2 qI 
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The detailed form of the coefficients of equation (|3J) is 

hl2l[l2 



( h I'l Ai 

h 1'2 A2 I . (57) 
V A'l A'2 L 



The range of the summation variables are determined by those of the coefficient H^^^^^i and those of 
the 9j symbol. 
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